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Abstract 

In this paper we prove the existence and local uniqueness of stationary states for 
the nonlinear Dirac equation 

3 

i j-'djip — mip + F(jlj'ip)'ip = 

3=0 

where m > and F{s) — \s\^ for 1 < 9 < 2. More precisely we show that there exists 
£0 > such that for cu £ {m — eo, m), there exists a solution ip{t, x) = e~^'^*<j)uj{x), xq = 
t,x = {xi, X2, X3), and the mapping from ui to (j>t^ is continuous. We prove this result by 
relating the stationary solutions to the ground states of nonlinear Schrodinger equations. 

1 Introduction 

A class of nonlinear Dirac equations for elementary spin-i particles (such as electrons) is 
of the form 

3 

i^-f^djip -71111; + F{-tpTp)'tp = 0. (1.1) 
3=0 

Here : M ^ M models the nonlinear interaction. ip : ^ is a four-component 
wavefunction, and m is a positive number, dj = d/dxj, and 7-' are the 4x4 Dirac matrices: 



h \ k f a 



where c'^ are Pauli matrices: 



^=10' " =[^ ' ^ = -1 
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We define 

2 4 

i=l i=3 

where (•, •) is the Hermitian inner product in C^. 

Throughout this paper we are interested in the case 

F{s) = \sf, O<0<oo. (1.2) 

The local and global existence problems for nonlinearity as above have been considered in 
[HE]. For us, we seek standing waves (or stationary states, or localized solutions of (jl.ip ) 
of the form 

V'(xo,x) = e-''^^<f>{x) 
where xq = t, x = (xi, X2, X3). It follows that </> : — > solves the equation 

3 

i Y l^djcj) -m(t) + uj-i^cj) + F{^(l))(j) = 0. (1.3) 

Different functions F have been used to model various types of self couplings. Stationary 
states of the nonlinear Dirac field with the scalar fourth order self coupling (corresponding 
to F{s) = s ) were first considered by Soler [iTj proposing them as a model of extended 
fermions. Subsequently, existence of stationary states under certain hypotheses on F was 
studied by Cazenave and Vazquez [2], Merle |5j and Balabane [T], where by shooting method 
they established the existence of infinitely many localized solutions for every < u; < m. 
Esteban and Sere in [3], by a variational method, proved the existence of an infinity of 
solutions in a more general case for nonlinearity 

m = Imr + ^i'A7Vr^)</', 7' = iViV 

for < ai,a2 < ^- Vazquez [15] prove the existence of localized solutions obtained as a 
Klein-Gordon limit for the nonlinear Dirac equation {F{s) = s). A summary of different 
models with numerical and theoretical developments is described by Ranada |10] . 

None of the approaches mentioned above yield a curve of solutions: the continuity of 
(j) with respect to uj, and the uniqueness of (j) was unknown. Our purpose is to give some 
positive answers to these open problems. These issues are important to study the stability 
of the standing waves, a question we will address in future work. 

Following [12\, we study solutions which are separable in spherical coordinates. 



(j){x) 
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where r = \x\,{9,^) are the angular parameters and f,g are radial functions. Equation 
(jl.3p is then reduced to a nonautonomous planar differential system in the r variable 

f' + lf = i\9'-fr-im-u^))a ^^^^^ 
g' = {\g^-fr-{m + u;))f. 

Ounaies in [H] studied the existence of solutions for equation (jl.Sp using a perturbation 
method. Let £ = m — u. By a rescaling argument, (II. 4p can be transformed into a perturbed 
system 

2 

r ' ' (1.5) 



v' + 2mn — e(l + \v 



2 £7.2 



If £ = 0, (jl.Sp can be related to the nonlinear Schrodinger equation 

At; , ,9/1 f' , , 

- — z;2% = 0, u=- — . 1.6 
Im 2m 

It is well known that for E (0,2), the first equation in ()1.6p admits a unique positive 
solution called the ground state Q{x) which is smooth, decreases monotonically as a func- 
tions of |x| and decays exponential at infinity (see [9], [13] and references therein). Let 
Uq = (Q)~2m^')' ^^^^ want to continue C/q to yield a branch of bound states with 
parameter e for (jl.Sp by contraction mapping theorem. 

Ounaies carried out this analysis for < < 1 and he claimed that the nonlinearities 
in p.5p are continuously differentiable. But with the restriction < < 1 we are unable 
to verify it. The term Iv"^ — eu'^\^ has a cancelation cone when v = zty^u. Along this cone, 
the first derivative of \v'^ — £u'^\^ is unbounded for < ^ < 1. But Ounaies' argument may 
go through for > 1, which gives us the motivation of the current research. However we 
can not work in the natural Sobolev space //-"^(M^, M^). Since ff-'^(M^) ^ L^(M^), we lose 
regularity. To overcome these difficulties, we want to consider equation (jl.Sp in the Sobolev 
space W^^'P(M3,M2),p > 2 and 6* > 1. 

To state the main result, we introduce the following notations. For any 1 <p < 00, Lr = 
Lr(M'^) denotes the Lebesgue space for radial functions on 'M? . Wr'^ = Wr'^{^^) denotes 
the Sobolev space for radial functions on M^. Let Xr = Wr'^ x Wr'^,Yr = Lr x Lr- Unless 
specified, the constant C is generic and may vary from line to line. In this paper, we assume 
that 771 = ^, since after a rescaling Tp{x) = (2m) ae xl'(2mx), equation (jl.ip becomes 
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We prove the following results: 



2 

j=0 
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Theorem 1.1 Let e = m — uj. For 1 < 9 < 2 there exists Eq = £o{0) > and a unique 
solution of ([HI) {f,g){e) G C((0, eo), W^r 1^^)) satisfying 

f{r) = e'4r{-Q\^r) + e2{V^r)) 
g{r) = e^{Q{^/er) + ei{^r)) 

with 

||ej||^yi,4 < Ce for some C{6) > 0,j = 1,2. 

Remark: The necessary condition < m must be satisfied in order to guarantee the 
existence of locahzed states for the nonfinear Dirac equation (see [15] , [6] ) . 

The solutions constructed in Theorem 1.1 have more regularity. In fact, they are classical 
solutions and have exponential decay at infinity. 

Theorem 1.2 There exists C(e) > 0,cr(e) > such that 

\ej{r)\ + \drej{r)\ < Ce"'^'^ j = 1,2. 

Moreover, the solutions {f,g) in Theorem 1.1 are classical solutions 

2<p<+oo 

Remark. From the physical view point, the nonlinear Dirac equation with F{s) = s (Soler 
model) is the most interesting. In fact, Theorem 1.1, Theorem 1.2 are both true for the 
Soler model. In fact, from (|1.5|) one can find out that (f^ — eu^) — = —ev? which is 
Lipschitz continuous. An adaption of the proofs of the above theorems will yield: 

Theorem 1.3 For the Soler model F{s) = s, there is a localized solution of equation (II. 3p 
satisfying Theorem 1.1 and Theorem 1.2. 

Next we proceed as follows. In section 2, we introduce several preliminary lemmas. In 
section 3, we give the proof of Theorem 1.1, Theorem 1.2. 

2 Preliminary lemmas 

We list several lemmas which will be used in Section 3. 

Lemma 2.1 Lei 5 : M ^ M 6e defined by g{t) = jip^t, 6* > 0, then 

g{a + 0-) - g{a) - {26 + l)|apV 
where Ci , C2 depends on and Ci = if < < ^. 



1 a 



\2e-i 



2F 



|2e-l^ 



-7 
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Proof. We may assume that a > in our proof. It is trivial if cr = 0. So we assume that 
(7 7^ 0. If a < 2\a\, then |a + cr| < 3\a\ and 

g{a + a)- g{a) - {29 + l)a^^(T 
< \g{a + a)\ + \g{a)\ + {2e + l)\a^'a\ 



If a > 2|cr|, then 
so that 

Taylor's theorem gives 



a + a> 2\a\ + cr > |a| > 0, 
g{a + a) = (a + af^+K 



g{a + a)- g{a) - {29 + l)a''a = -/(^a^ 

where ^ is between a + a and a. Since g"{(,) = 29{29 + if 26* - 1< 0, then 

1/(01 < C\af'-\ 

If 26* - 1 > 0, we have 

1/(01 < Cmeo,{{a + af'-\a^'-'} < C{\a\^'-' + lal^^"^). 
Hence we prove the lemma. □ 

Lemma 2.2 For any a,b G R,9 > 0, we have 

|a-6|^-|a|^ < Ci\af-^\b\ + C2\bf 

where Ci , C2 depends on 9 and Ci = 0if0<9<l. 

Proof. The proof is basically similar to that of the lemma as above. It is trivial if 6 = 0. 
So we may assume that 6/0 and a > 0. li a < 2\b\, then 

\a-bf-\af <C(|a|^ + |6|^) <C|6|^. 

On the other hand, if a > 2\b\, then |a — 6| > a — \b\ > \b\. So by using the mean value 
theorem 

\a-bf-\af =9\tf-^\b\ 
where t is between a — b and a. If — 1 > 0, then 

\tr'<c{\af-'+\br'), 
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hence 

|a-6|^-|a|^ <Ci\a\^-^\h\ + \b\^. 
If 6* - 1 < 0, then \t\^-'^ < C|6|^"\ so that we conclude 

|a-6|^-|a|^ <C|6|^ 

The proof is complete. □ 

Lemma 2.3 For any a, 6, c S M, if 1 < 9 < 2, then 

\a + b + cf -\a + bf - |o + cf + \af < C{\cf~^ + \bf~^)\b\, 
where C depends on 9. 

Remark. This inequality is symmetric about 6, c, so the right hand side can be equivalently 
replaced by C(|c|^~^ + |6|^~-'^)|c|. Without loss of generality, we assume that |6| > |c| in the 
following. 

Proof. For simplicity, let 

L = \a + b + cf -\a + b\^ -\a + + \a\^ . 
It is trivial for 9 = 1, since if \a\ > 5|6| then L = 0. If \a\ < 5|6|, 

\L\<C{\b\ + \c\). 

So next we consider > 1. If |a| < 5|6|, by triangle inequality and Lemma 12.21 we have 

1^1 <C{\a + c\'^-^ + \af-^ + \bf-^)\b\ 
<C{\cf'' + \bf-')\b\. 

If |a| > 5|6|, by using Taylor's theorem 

\L\=C\{a + tib + t2c)f-^\bc\. 

where ti,t2 G (0,1) and 

|(a + hb + t2c)\ > \a\ - 2\b\ - \c\ > |c|. 

So if 1 < < 2, we have 



\L\<C\cf-'\b\. 



The proof is complete. 



□ 
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Lemma 2.4 Let 2 < p < oo, / : — M 6e radial and bounded. Suppose fr + € 
^foc f/ e i^L- ^//^ + f / e i/ien ieLP and 

\\-\\Lv<C\\drf+-f\\Lv 

r r 
Proof. We begin with p = oo. Using integration by parts 

V(r)= f{dpf + -f)p^dp. (2.1) 



^2. 



Hence ^ 

|rVl < Wfr + -/IIl- r s'^ds = U\fr + -/IIloo 

r Jo 3 r 

which gives 

II^IIloc <C||5,/ + ^/||loo. 

Next let us consider p = 2. Let < ri < r2 < oo. Denote D = {x e R^,0 < ri < |x| < r2} 
and 

Jri r r 

By Holder inequality, 



/ 2 

-'"<C'l|-||L2(I3)ll/r + -/||L2(D)- 



On the other hand, we have 



r = l\\^\\hiD)+^\r2f{r2)-r^f{n)). 
Since r2p{r2) > we have 

ll^lli^p) < C (\\^\L^^D)\\Ur + lf)\\L^^D)+nf\n) 



Let r2 — ^ oo, ri — > 0, we obtain 



<c||9, + ^/IU2. 



The intermediate case 2 < ^? < oo is a direct result of interpolation . □ 
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3 Proof of the main theorems 

Similar to [8], we use a rescaling argument to transform (|1.4p into a perturbed system. Let 
e = m — uo (remember m = ^). The first step is to introduce the new variables 

f (r) = u{\/er) , g{r) = e^v{\/er) 
where {f,g) are the solutions of (jl.4p . Then {u,v) solve 



u' + -U- \v\^% + v- (\v^ - en¥ - \v\^'^)v = 0, 
r 

v' + u-e{l + \v'^ - eu'^f)u = 0. 
Our goal is to solve (|3.ip near e = 0. If e = 0, (|3.ip becomes 



(3.1) 



u' + -U- Ivl'^^v + v = 

r (3.2) 

v' + u = 0. 

This yields the elliptic equation 

-Av + v = \v\^%, u = -v' (3.3) 

It is well known that for < ^ < 2, there exists a unique positive radial solution Q{x) = 
Q{\x\) of the first equation in (13. 3p which is smooth and exponentially decaying. This 
solution called a nonlinear ground state. Therefore Uq = {—Q', Q) is the unique solution to 
()3.3p under the condition that v is real and positive. We want to ensure that the ground 
state solutions Uq can be continued to yield a branch of solutions of (|3.ip . 
Let 

v{r) = Q{r) + ei(r), u{r) = -Q\r) + e2(r). 
Substitution into (13. ip gives rise to 



e'2{r) + -e2(r) + ei - {26 + l)Q'^ei = Ki{e, 61,63) 

r (3.4) 

ei(0 + e2(r) = -R'2(e, 61,62) 



where 



i^i(e,6i, 62) = IQ + 6i|2^(Q + ei) - {26 + l)Q2eei - Q^e+i 

+ (|^2_^^2|e_^2e)^ 

i^2(e, 61,62) = e(l + 1^^ - ev?\^)u. 
Define L the first order linear differential operator L : Xr Yr by 

ex \ _ (\-{26 + l)Q2e 5^ + I\ / 



' 62 y V 1 y V 62 
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Then we aim to solve the equation 

Le = K{e,e) (3.5) 

where e = {ei,e2)^,K{e,e) = {Ki,K2)'^{e,e). Let / = {0,a),a > 0. We say e{e) is a weak 
X^-solution to equation ()3.5p if e satisfies 

e = L-^K{e,e) (3.6) 

for a.e. e £ I. L is indeed invertible as we learn from the following lemma. 

Lemma 3.1 Let < < 2, the linear differential operator 

\ dr 1 

is an isomorphism from Xf onto Yr for 2 < p < oo. 

Proof. First we prove that L is one to one. Suppose that there exist radial functions 
ei, 62 G Wr'^ such that 

1=0. 



62 

Then 

- Arei + ei- {29 + l)Q^^ei = 0, e2 = -e[. (3.7) 

It is well known (see, eg. |14j ) that ei = is the unique solution in H^. 

Next we prove that L is onto. Indeed L is a sum of an isomorphism and a relatively 
compact perturbation: 

M is relatively compact because of the exponentially decay of the ground state at infinity. So 
we only need to prove that L is an isomorphism from Xr to Yr, i.e. for any (0i, ^2) G LrXLr, 
there exist (61,62) G Wr'^ x Wr'^ such that 

I I \ _ { 'Pi 

It is equivalent to solve 



62 / \ 92 



ei + {dr + ^)62 = 01 

a^ei + 62 = 02 

and show that 61, 62 G w}'^ . By eliminating 62 we know that 61 satisfies 

(-A, + 1)61 = 01 - (5. + -)02. (3.9) 

r 
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Define G{x) = (47r)- V|~^e~i''i. I^Mj has the solution 

ei = G{x) * f 01 - {dr + -)(t>2 
\ r 

= G{x) *</>!+ drG{x) * (t)2- 

Here we have used the property of convolution and the fact {dr + ^Yfir) = —drf{r) in 
By Young's inequality and G,drG £ L^(E^), we have 

lleiliLP < \\G\\L4MLr + \\drGh4<P2\\L^^ 

which implies 

ei G LP. 

Similarly 62 satisfies 



2 

-Ar + 1 + — )e2 = (p2- drCpl- 



Let = j^G{x), then 



62 = -fr(x) *(j)2- H{x) * {dr<l)l) 

= Hix) * (/>2 + idrH + * 01 G LP 

r 



since if, (9^ + G 

To improve the regularities of ei, 62, we go back to ()3.8p . Since 

drei = 01 - 62 G L^, 

we have 61 G Tl^r'P. Regarding the regularity of 62, we know that 

2 

{dr + -)62 = 01 - 61 G LP. 
r 



By Lemma [27 

||a.62|U. <C(||^||iP + ||(a, + ^)62||,..) 

<C||(5, + ^)62|Lp = C||0i-ei||^P. 
Hence we have 62 G VTr □ 



Now we are ready to construct solutions of (13. 6p by using the contraction mapping 
theorem. 

Proof of Theorem 1.1. To prove Theorem 1.1, we prove there exists ei > such that 
for every < e < Si, there is a unique solution to equation p.6p 

6 = L-^K{e,e) 
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in a small ball in X^. First we must ensure that K{e,e) is well defined in Yr if e G Xr- 
Recall that 

i^i(e,ei, ea) = \Q + ei\^\Q + ei) - {29 + l)Q^^ei - Q^^+^ 

+ {\v^ - eu'^f - v^^)v 
K2{e, 61,62) = e{l + Iv"^ - eu^f)u. 

Let us consider Ki, the estimate for K2 is similar. Since 

\K,{e,6)\<C,,e{\v\''^' + W\''^') 

where C^^e is a real constant depending on e,6, it suffices to show that dt'p^"'"^ + S 
LP. By Sobolev's embedding 

for any g if p > 3. We choose p = 4 in the following. The same argument is available for 
K2. From Lemma |3.H we know that L'^K £ X? . 
Fix (5, to be chosen later. Consider the set 

^ = {6&XU\\6\\x^<5}, 

and suppose e G il. We know that 

\\L-^K{e,6)\\xf < Ci\\K,{s,6)h. + ||i^2(^, 6)1^4). 
Let Ki{e, e) = Kf{e, e) + isTf (e, e) where 

K^e, 6) = \Q + ei|2^(Q + ei) - (29 + l)Q^'6i - Q^^+i 

and 

Klie, 6) = (|(Q + 6if - ei-Q' + ea) Y " l(Q + eOP'') (Q + ei). 

Thus 

l|i^l||L4<ll^fllL4 + ||i^rilL4. 

For ||iCi||i4, let a = v'^ = {Q + 6if ,b = ev? = e{-Q' + 62)^ in LemmaE^l then 
\\KI\\l, < Cee \Q + eip^-^l - Q' + + | - Q' + e2p^|Q + ei| 



Li 



< CeeiWQWZP + ll^ll^r) 
<Qe(||Q||^tl+5)< 5/4 

if 5 < 1 and e is small enough such that 



Cee{\\QfZP+S)<^-. 
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For ||JC"||j;^4, let a = Q{r),a = ei in Lemma [2?T| then 

<C,(||eif^M + ||ei||j;tl) 

< Ce{d^ + < 2Cg6^ < 5/4 

if S < A similar argument can be applied to K2 (with similar condition on e,6) to 
obtain that 

\\K2{e,e)h.<i 



Hence we obtain 



Next we want to show that for any e, f £ ^l, and 5,e as above, 

\\L-\K{e,e) - KieJMy. <l\\e- f\\j,r, 

i.e. L~^K is a contraction mapping. We have 

\K{e, e) - K{e, f)\ < {KUe) - KUDl + \Kt{e) - Kt{f)\ + \K,{e) - K^if)]. 
We compute the r.h.s. term by term. After rewriting i^f (e) — K^{f), 

iK^e) - Knf)\ <\\Q + ei\^\Q + ei) - \Q + M^'iQ + h) - {29 + l)\Q + h\^\ei - h) 
+ {26 + 1) I IQ + /i|2^(ei - /i) - Q^%ei -f^)\=D^ + D^. 

For D", let a = Q + /i, o" = ei — /i and by use of Lemma [2m then 



i?r<c(|g + /i|2^-i + |ei-/i| 



By Sobolev embedding and Holder inequality, we have 

11^1 < Cg{\\e, - + ||ei - hW^P 



<Ce{5 + 5'')\\e^-h 



< il|ei-/illwi^4 



if (5 < j^Q-- Using Lemma [221 we find 



1^2"! < Ce{Q^'-^ + \2Qh + /i |'"')|2g/i + /fllei - /i| 



Hence 



\\mLl<CB5\\e^-h 



,v < -||ei - /i 
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Then let us study Kf{e) - KK f) : 

Klie) - Klif) = (\iQ + ei)2 - e(-Q' + es)^ - HQ + eiff) (ei - /: 



+ (|(Q + ei)2 - e{-Q' + 62)^ - HQ + eiff){Q + /i) 
- (I(Q + /i)' - e{-Q' + /2)Y - l(Q + /i)T)(Q + /i)- 
Notice that the first hue in the r.h.s. is easy to estimate since 

){Q + eif - e{-Q' + e^ff - \{Q + e^ff) {e, - h) 



< Cee 
1 



{\Q + eiP'-'l -Q' + + \-Q' + 62^(6, - h) 
< gllei -/illi^^.p 

for e sufficiently small. For the second and the third line, let us define 

Eie, f) = (|(Q + ei)2 - e{-Q' + 62)^ - \{Q + ei)Y)(Q + /i) 
-{\{Q + h? - e{-Q' + /2)T -\{Q + fi?f){Q + h)- 

We discuss the contractive property for two different situations > 1 and 9 = 1 separately. 
For e>l,we use Lemma O Set a = {Q + f^f,h= {Q + eif - (Q + fif ,c = -e{Q + f2f 
(notice that 6, c can be taken sufficiently small), and rewrite E{e,f) to get 



\E{eJ)\ < 
+ 



\a + b + c[ 



\a + bf 



a + cr + a' 



\{Q + eif - e{-Q' + e2?\' -\{Q + ex? - e{-Q' + f2? 



<CB{\b\'-' + \cf'^] 



(3.10) 



+ 



\{Q + ei)2 - ei-Q' + 62)^ - |(Q + e,f - e{-Q' + f2? 
where for the last line, we applied Lemma |2. 21 We obtain 

\\E{e, f)\\Lf. < Ce{e'-' + 5'-')\\e - f\\^y <\\\e - /||^,4 
for e, 5 sufficiently small. Hence we have for 1 < 9 < 2, 

\\K,{e)-K,{f)\\Lf<l\\e-f\yiA. 

Next we prove that E{e, f) is contractive for 9 = 1 directly. Lemma 12.31 can not be used 
since |6|^-^ = |c|^-^ = 1. In (l3lTO . if |o| > max{5|&|, 5|c|}, then 



|a + 6 + c| — |a + 6| — |a + c| + |a| = 0. 
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Thus 

\\E{e, f)\\Lf < Cee\\e - /||^^m <\\\e - f\\^y. 
Hence we only need to consider E{e,f) if \a\ is small, i.e. if \a\ < 5max{5|6|, 5|c|}, 

\a + b + c\-\a + b\-\a + c\ + \a\ < C(|6| + |c|). 
Simply assume that |c| < we have 

\\Eie, < Ce6'/^e - f\\w^A <\\\e- fWw^A. 

Therefore if ^ = 1, 

\\K,{s,e) - K,ie,f)\\L. <^\\e - fW^y. 
Similarly, we can prove that 

1 

4 



K2ie)-K2if)\\y.<-\\e-f\\^.„. 



Note we can satisfy all the condition above by choosing 6 = Cge and taking e sufficiently 
small. Then the contraction mapping theorem implies L~^K has a unique fixed point 
e(e) G O which is a weak solution of equation (j3.5p . The continuity w.r.t. e follows from 
the continuity w.r.t e of the map L^^K and its contractibility. This completes the proof of 
Theorem 1.1. □ 

Let us see why a solution of equation (j3.5p which is in has more regularity. This is 
done by using a standard bootstrap argument and the following standard lemma: 

Lemma 3.2 Let F : C ^ C satisfy F{0) = 0, and assume that there exists a > such that 
\F{v)-F{u)\<C{\v\'' + \un\v-u\ for all u,veC. 

Let 

1 a 1 

- = h -, I < p,q,r < oo. 

r p q 

It follows that ifu€ LP, Vn G L'^, then VF{u) e L'' and 

\\VF{u)\\Lr < C7||n||2p||Vn||L,. 
Proof of Theorem 1.2. First we can prove that 

ei,e2G fl W^'P. 

4<p<oo 



Recall that ei , 62 satisfy 



e2 ~ \K2 
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and L is an isomorphism from X"^ YJ. We know that 
and 

<e(|^| + |^|2m + |^|2m) 

Since 61,62 S Wr'^^iM?) ^ then Ki x K2 e L? x for p > 4. By LemmaEH we 

have 

(ei,e2)e fl W^'^xW^'P. 

4<p<+oo 

Next from Lemma 12.11 Lemma 12.21 and Lemma 12.31 

|i^(6i, 62) - /2)| < C(|Q + Qf' + \e,\'' + |62|2^ + \h\'' + 1/21''^) 

(|ei-/i| + |e2-/2|). 



So by Lemma 



This gives that 



and 



4<p<oo 

(ei,e2)G fl W^'PxW^'P 

4<p<oo 



Going back to equation (j3.5|) . we know that 61,62 G Wp^ C C^. So (/,c/) are classical 
solutions. 

Moreover we show that 61,62 have exponential decay at infinity. We know 61,62 are 
classical solutions and |6i|, I62I < Ce by Sobolev's embedding theorem. Taking derivatives 
in (j3.2p and after tedious computations we find 

61" - 61 = 5i(r)6i + (52(r)6'i + 53(r)Q for r large (3U) 
62" - 62 = o-i(r)62 + (72{r)e'2 + cr3{r)Q for r large 

where ai,Si G W'^'^ and |cri|,|(5i| < Ce{i = 1,2,3) for r large. 

We conclude that there exist constants ro, C(e) positive such that 

|ei(r)| + |62(r)| < Ce'"'' for r > vq. (3.12) 

We prove it by an application of the maximum principle. Without loss of generality, suppose 
ei(^o) = 2e {tq is sufficiently large). Let 

h{r) = e-^(''-'^«) + /36'^('-^o) 
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where /? > is arbitrary and < z/ < 1 is to be determined later. If g = ei — h, then g 
satisfies 

/ = (1 + 6i)g + 62g' + (1 - i/^ + 5i)h + 62h' + 63Q 
Since /i' = i/(-e-''('"-'^o) + /Je'^('^-^'o)) < jy/i and Q < h, then 

/ > (1 + Si)9 + -^a^' + (1 - + <5i + |53|)/i (3.13) 
with ^(ro) = ei(ro) — (1 + /?) < 0, ^(cxd) < 0. Thus we claim that 

g{r) < for r > tq, 

if v is small enough such that 

l-iy^ + 6i + \53\ >0. 

If the claim is not true, then g(r) obtains maximum at r = ri and g{ri) > 0. Thus 
g"{ri) < 0,g'{ri) = 0. But this contradicts with equation (|3.13|) since the right hand side of 
(j3.13p is positive evaluated at r = ri. Therefore the claim is true if v < \/l — Ce and then 

ei(r) < h{r) if r is large enough. 

Then similarly we can show that 

> —h{r) if r is large enough. 

Thus 

|ei(r)| < h{r) = e'^^''-'"°'^ + /3e'^('^-'^«). 

Letting /5 ^ , we have 

|ei(r)| < Ce-''''. 

for r large enough. The exponential decay estimate for 62 can be obtained in a similar way. 
Once we have (]3.12p . it is obvious that \drej{r)\ < Ce"*^^ and Cj E H^. This completes the 
proof of Theorem 1.2. □ 

Remark. For < ^ < 1, our method does not work since Lemma 12.31 is not valid. Let us 
consider a special example. Suppose 62 = /2 = 0, then 

E{eufi) = {\{Q + eif - e{Q'f\' -\Q + e^^' - |(Q + hf - e{Q'ff + IQ + h?'){Q + /i) 
We want to know whether or not the following inequality is true 

l^(ei(r),/i(r))| < ^|ei(r)-/i(r)|, r G (0, 00) (3.14) 

if e small enough. Letting tq large enough and s = > to be determined later, we 

assume that 

Qiro) + ei{ro) = ^\Q'{ro)\il + s), 
Q{ro) + fi{ro) = V^\Q'{ro)\. 
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Then under this ansatz, 

\E{ei{ro),fi{ro))\ = [{s' + 2s)' - ((1 + sf' - l)]h''+' = g{s)h^'+\ 
|ei(ro) - /i(r-o)| = sh 

where h = -v/e|Q'(^o)l- Then 

\E{e,{To)Ji{ro))\ = ^/i2^|ei(ro) - /i(ro)|. 
We claim that if a > j-^, then 

^/i2^>-, as e^O. 
s 2 

g{s) > Cs' 

\{l + sf' -l\ < C7(s + s^^) « C/. 



In fact, we have 



since 



and 



So 

9{s) 



s - ' 2' 

since 9 + a{9 — 1) < 0. The claim is proved and consequently, ()3.14p does not hold for every 
r G (0, oo). 
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